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Q \ Abstract 

, We investigated the relation between the behavior of gravitational wave at late time 

and the limit structure of future null infinity tangent which will determine the topol- 
ogy of the event horizon far in the future. In the present article, we mainly consider 
£N) ■ a spacetime with two black holes. Although in most of cases, the black holes co- 

^ ' alesce and its event horizon is topologically a single sphere far in the future, there 

are several possibilities that the black holes never coalesce and such exact solutions 
as examples. In our formulation, the tangent vector of future null infinity is, under 
conformal embedding, related to the number of black holes far in the future through 
the Poincare-Hopf 's theorem. Under the conformal embedding, the topology of event 
horizon far in the future will be affected by the geometrical structure of the future 
f-— ^ null infinity. In this article, we related the behavior of Weyl curvature to this limit 

behavior of the generator vector of the future null infinity. We show if Weyl curvature 
decays sufficiently slowly at late time in the neighborhood of future null infinity, two 
black holes never coalesce. 
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1 Introduction 

> ■ 

We are often interested in the final state of a black hole spacetime after gravitational collapse. When 

S— i ' 

matter collapses, it passes inside an event horizon and energy will be emitted to future null infinity in 
the form of gravitational waves. Since the amount of this energy is limited, we may expect that the 
spacetime will approach to a stationary state. If a black hole spacetime becomes stationary far in the 
future to an observer outside the black hole, the final state of the black holes is considerably restricted. 
As used in reference |10| . the term far in the future means sufficiently at late time to an observer which 
remains outside a black hole. 

So far the final state of the black hole spacetime has been studied by many authors and these studies 
are known as uniqueness theorems of black holes. Israel showed that the only static and topologically 
spherical black hole is the Schwarzschild solution pQ or the Reissner- Nordstrom solution 2 . In stationary 
axisymmetric situations, the uniqueness theorems of a vacuum black hole and a charged black hole were 
shown by Carter Robinson an d Mazur 

All they are based on the topological notion, which is also well studied for stationary black holes[H]|7][S]. 
Hawking proved that the spatial topology of the smooth stationary horizon must be spherical 6;. More- 
over, the work of Chrusciel and Wald 7 showed that in a stationary spacetime, the spatial topologies of 
the connected components of black holes are only spheres under null energy condition. 
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Hence, now it is known that, under some reasonable conditions such as asymptotic flatness and 
the weak energy condition, each component of the black hole region is topologically trivial, i.e., simply 
connected^. On the other hand, there were numerical simulations which suggest non-trivial topologies 
of the horizons[D]. There have been some confusion, but it is now well understood that even though 
the black hole region in the spacetime is simply connected, there are many possible topologies of spatial 
sections. In particular, one of the authors^H] showed how topology of the spatial sections of a black hole 
is related to the endpoint set, or similarly, the crease set, of the event horizon. 

Thus, we understand that the topology of the event horizon can be realized by the formation stage 
with the crease set and the final stationary stage, independently. If the black hole spacetime settles 
down to a stationary state, we expect the spatial topology of the event horizon will consist of some 
connected components with spherical topology. To realize these aspects more, a new approach is required 
to investigate the topology of an event horizon far in the future. Especially the number of connected 
components of the black hole is most important. 

There is not a simple method to know topological structure of an event horizon far in the future, that 
is, to determine how many black holes the final state of the black hole spacetime will consist of. Then 
we consider some examples, they are Majumdar-Papapetrou solution|14|. C-metric|15| (which recently 
becomes important in studies of higher dimensional black holes since the black ring solution j 1 1 j with 
vanishing angular momentum corresponds to a higher dimensional C-metric) and so on. The spatial 
topology of the event horizon of such a solution is obviously not a single sphere. 

Since in an asymptotically flat spacetime;! 7. , (to be more exact, in a strongly asymptotically pre- 
dictable spacetime) an event horizon is defined as the boundary of the causal past of future null infinity 
,f + , we expect the topology of an event horizon will be strongly related to the topological structure of 
future null infinity. Because the topology of future null infinity in an asymptotically flat spacetime is 
S 2 x R, we might conclude far in the future the spatial topology of an event horizon will always be a 
single sphere. Nevertheless it is well known that there are several above mentioned stationary solutions 
including black holes separated from each other. 

Therefore we have to study the structure of the neighborhood of in detail, especially in the limit 
of the future direction. The geometric study of the asymptotic region was successful in the formulation 
of large distance expansion (the expansion of a conformal factor O of conformal embedding) like the 
investigation of the peeling property j!2j. Then we evolve the large distance expansion of geometry into 
the future direction along the tangent of the future null infinity n a = g ah Vi,f2. 

To analyse that, a key aspect is that there remains a long tail of Weyl curvature ^4 in the C- 
metric spacetime, which is discovered by Tomimatsu 13 . In the present article, we develop pseudo- 
Newman Penrose basis vectors and clarify the relation between the Weyl curvature component of outgoing 
gravitational wave and the congruence of the integral curve of the tangent vector n a . We will find 
these aspects rule the topology of an event horizon in the limit of those integral curves. 

To determine the gauge condition of the conformal embedding for studying this limit structure of the 
tangent field n a , we should be careful since such a limit will be realized only when conformal embedding 
makes unphysical manifold compact in this direction. We show from the discussion in such a unphysical 
manifold that two black holes are separated forever for an asymptotic observer outside black holes, if its 
Weyl curvature in the neighborhood of future null infinity decays slowly at late time. 
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Since we discuss black hole spacetimes, only asymptotically flat spacetime is considered. So, in section 
121 we will give the notion of an asymptotically flat spacetime and the gauge condition of conformal 
embedding. The congruence of the integral curves of tangent field n a — g afc V fc f2 is investigated in the 
third section. Fourth section provides the index theorem and using it, study the relation between the 
topology of an event horizon and the structure of that congruence. In section we will apply this 
discussion to C-metric as an example, to illustrate its detail concretely. The final section is devoted to 
summary and discussions. Throughout this article, we use the abstract index notation as the component 
notation of tensors and it is denoted by Latin indices a, b, ■ ■ ■ . 

2 Asymptotic flatness and gauge choicing of conformal com- 
pact ificat ion 

In this article, we use the term far in the future to mean at sufficiently late time to an observer outside 
black holes in stead of referring the time after a spacetime settles down to a stationary phase. Then, our 
main question is "How can we know the topology of an event horizon far in the future ?" 

To discuss the topology of an event horizon, we pay attention to the relation with future null infinitiy 
since the event horizon is defined as the boundary of the past set of future null infinity H + = d( J~ (J? + ))D 
M. Because of its non-compactness (^ + ~ S 2 x R)|17| there is an ambiguity of the limit structure. 

Here our tactics is to compactify the future null infinity and the integral curves of n a tangent to 
+ in their direction to discuss the limit of the integral curves. Therefore inextendible endpoints of the 
integral curve of n a will be attached after gauge choice of the conformal embedding. It will be discussed 
that a geometric condition to allow that a set of endpoints with non-trivial topology is attached, by the 
Poincare-Hopf theorem. From this we discuss the topology of an event horizon far in the future. 

First of all, we mention asymptotic flatness. Asymptotic flatness is essential in our discussion. We 
have to state the definition of asymptotically flat spacetime before explaining the analytic method to 
examine the limit structure of J r+ . The definition mentioned here is based on the reference |17j . 

Definition 1 (Asymptotic Flatness). A spacetime (M, g a b) is said to be asymptotically flat at null 
infinity if there exists an unphysical spacetime (M, g a b) with g a b = 2 g a b smooth everywhere except at the 
point i° where it is C >0 and conformal isometry ip ; M — > ip(M) £ M with conformal factor 51 satisfying 
the below conditions. 

1. J + (i°) U J~(i ) = M — M . Thus, i is spacelike related to all points in M and the boundary of M 
consists ofi°, J? + = j + (i°) - i° and J~ = J~(iP) - i° . 

2. There exists an open neighborhood V of M = i° U U such that the spacetime {V,g a b) is 
strongly causal. 

3. Q, can be extended to a function on all M which is C°° everywhere and C >0 at i . 

4. On J + and we have = and V a ^ 0. 

5. The map of null direction at i° into the space of integral curves of n a = g ab Wb^ on J? + and is 
diffeomorphism. 
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6. For smooth function on M — i° with > on M U U which satisfies V a (u> 4 n a ) — on 
J^ + , the vector field <jj~ x n a is complete on .JP^ and ^~ . 

7. In the neighborhood of J? + and J?~ physical Ricci tensor behaviors as R ab = 0(fl 2 ). 

Let us note that there is gauge freedom in the choice of an unphysical spacetime (M, g ab ) with an 
asymptotically flat physical spacetime (M,g ab ). This gauge freedom is important in our discussion. 

// the spacetime (M, g ab ) is an unphysical spacetime satisfying the above definition with conformal 
factor 0, so is (M, u> 2 g a b) with conformal factor 00 for the any function to which is smooth everywhere 
except at i and positive everywhere. 

The most common and useful gauge condition will be an affine gauge. In this gauge, the tangent 
n a = g ab Vb£l of J? + generators is affinely parametrized. See the equation (11.1.21) in reference |17|. 
We should note that in this gauge, the null generators n a (so, the future null infinity generated out of 
them) are not compactificd in the future direction, since the integral curves of n a have infinite ranges 
of affine parameters in the future direction. Nevertheless in the present article, we consider another 
gauge condition such that the generator of can be compactified. Our purpose is to investigate the 
congruence of the n°-integral curves around the future null infinity in the limit of their parameter. To 
study this, it will be valid to compactify the future null infinity and the integral curves in their direction. 
The compactification will be realized by attaching any endpoints of the n a -integral curves, as + has 
been attached to M to study its asymptotic structure. Therefore u> vanishes at the set £ of attached 
endpoints. Taking a sufficiently small neighborhood of £, lo is considered to be monotonically decreasing 
in the direction of n°-integral curves. We call this gauge a compact gauge. 

Here it should be commented that according to the gauge condition, the integral curves go to 
i + or singularity. 

3 Limit structure of future null infinity 

We investigate the geometrical structure around future null infinity in a limit of the integral curves of 
n a = g ah VbVt. We define the limit as the limit approaching to the attached endpoints £. We call this 
limit rt a -limit. 

3.1 Physical Ricci Tensor and Unphysical Ricci Tensor 

The relation between the physical Ricci tensor and the unphysical Ricci tensor gives us a great deal 
of information. The physical Ricci tensor R a b is related to the conformally translated unphysical Ricci 
tensor R a b by 

R ab = R ab + 20- 1 V a 'V b 

+ga b g cd (0- 1 V c V d - 30~ 2 V c 0V rf 0). (l) 

By multiplying 2 and taking the limit — > we find that a vector n a = g ab Wb0 must be extended 
smoothly to J? + and be null at J !+ because the first term vanishes from conditions 4 and 7 of the 
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definition^ the second vanishes from conditions 3 and 4, and the fourth vanishes from the condition 3,4 
and g a b being smooth at J + . 

n a is null at ^ + but is not null off of J !+ and we put the shift from null vector as the following 

g ab n a n b = k (1) + k (2) 2 + ■ ■ ■ , (2) 

where ku\ (i = 1, 2, • ■ • ) are the functions independent of 0. The integral curves of n a on J> + is null 
geodesic generators of J? + , but the integral curves of n a is not geodesies off of J? + . Let us put the shift 
from the geodesic as 

n a V a n b -p (o) n b = O(0). (3) 
In particular, we put the contraction of the above equation with l a as the following: 

l b n a V a n b = -;p (0) + p (1) + ■■■ , (4) 

where we determine the normalization as l a n a = —1 (see the below equation l|10|l ). 
From and Q , we obtain the relation between k^ and pu\ 

1 3 

P(o) = 2 fc (!)' p W = " fc ( 2 )' p ( 2 ) = ~2 fc(3) ' " ' ' ^ 

where we used the relation of torsion free ; 

V a nf, = V b n a . (6) 
3.2 Pseudo-Newman-Penrose Formalism 

Newman- Penrose formalism JH] is the choice of the null basis which consists of a pair of real null vectors, 
l a , n a and a pair of complex conjugate null vectors m a , m a . Now we consider n a = g ab Vb0 as one of real 
null vectors in Newman-Penrose formalism. As we have mentioned before, n a is null on ^ + but is not 
null away off J? + . Therefore, strictly speaking, this is not Newman-Penrose basis vectors. Nevertheless 
most of Newmann Penrose framework are not affected by this difference. 
So the following orthogonality conditions are required 

l a m a = l a fh a — n a m a — n a fh a = (7) 

in addition to the conditions that vectors are null, 

H a = rn a m a = fh a m a = 0, (8) 
and the fact that n a = g ab Vb0 is null only on J^+, 

n a n a = k (1) + k (2) 2 + ■ ■ ■ . (9) 
We impose on the basis vectors the further normalization conditions, 

l a n a = -1, m a fh a = 1, (10) 

which is consistent with the definition of the vector n a and l a . Thus, the metric can be represented by 

g a b = -(fc(i)0 + fc(2)0 2 H )hh - 2l(a n b) + 2m (a m b) (11) 
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in this pseudo-Newman-Penrose basis vectors. Note that the first extra term of the above equation is 
caused by the equation |(?JJ|. 

In order to investigate the congruence of the n a -integral curves around J?^ ', we have only to consider 
the values of spin connections and curvatures only in the neighborhood of . Now we expand the each 
basis component of the equation Q in the powers of conformal factor as x = x^ +X( 1 ' j Q + X(2)Q 2 + .... 
The (n, I) component of the equation JQ) gives 

O(0 2 ) = Qfc(i)-3 P(O) +2 A1(O) ^0- 1 

+ fe? + 3p(x) + 2%) + 2/z (1) + k {1) qj + 0(0) (12) 

(13) 

where we used the condition 7 of the above definition ^ an d the equations 121,0 and fi = —m a rh b V a nb 
is the spin connection, which means the expansion of the integral curves of n a . q is the function satisfying 
l a VJ b = ql b . 

So we obtain from the coefficients of the first and the second term 

5 

-fc(i) - 3p( 0) + 2/i( ) = (14) 

and 

R-nl + 3 P(D + 2fc (2) + 2 M(i) + fc(i)5 - 0. (15) 

Similarly, from other components we can obtain the relations between spin connections and Ricci 
tensors. 



From the (1,1) component, we obtain 



From the (m, m) component, we obtain 



g = 0, (16) 



4 0) =0. (17) 



A(o) — A( ) — 0, (18) 



A(i) = RtL (19) 

where A = —rh a fh b V ' a rn, is the spin connection which means the shear of the integral curves of vector 
field n a . Eq^| means (ft = 0) is diffeomophic to S 2 x R. The (n,m) component gives 

f(o) = ^(o) = 0, (20) 
2i/ (1) = RW, (21) 

where v is defined by v = n a "fh b 'S/ a ni ) . 
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3.3 Change of Conformal Factor 

Now it should be recalled that there is gauge freedom in the conformal transformation considered above. 
Under the transformation 0-»0' = 00, g ao — > g' ab — 2 g a b, we have 

n a n' a = 0- 1 n a + 00~ 2 g ab W b 0. (22) 
In particular, on .y + the above transformation becomes 

^a la ^ — la /oo\ 

n — > n =0 n . (Z6) 
Under this transformation, the expansion of the null geodesic generators is transformed as 

Ho) "> Ko) = ^f-~ ^ na ^a0 (24) 

on . Therefore, note that we can choose /i(o) as any function independent of 0. Given A*/ )> smce the 
equation l(24|l is merely an ordinary differential equation, there always exists satisfying the equation 
(|24ll . From the equations (|5l) (|14l) . we find the relation 

M(o) = ~P(o) = -| fc (i)- ( 25 ) 

So from the equation |@J, we see that in the gauge satisfying ^( ) = 0; the nuu geodesies on y + are 
affincly parameterized. In this article, we call this gauge affine gauge. When we discuss the limit structure 
around , we should not choose the affne gauge since in the affine gauge, we cannot compactify J'^ 
or n a -integral curve and is not convenient to study the congruence of n a -integral curves in the n a -limit 
on neighborhood of J^ + . Hence, we choose the compact gauge such that J> + can be compactified by 
attaching a set £ of future endpoints of n a , that is, — > in n a -limit. 

From 1)18(1 . the leading order for the shear A, vanishes, which does not depend on the choice of the 
gauge; 

*(0) = A (o) = (26) 

holds. Then also in this gauge, spatial section of J !+ (fl = 0) is always S 2 . J'^ is topologically 
S 2 x (R U {+oo}) after future endpoints of its generators are attached. Under the transformation — > 
0' = 00, the first order is transformed as 

A 1 2 
A (1) -> A' (1) = -Si - —m a fh b V a V b - -(m a V Q0 ) 2 . (27) 

From now on, ' represents the quantity after transforming from an affine gauge into a compact gauge. 

3.4 Congruence of n a -curves and Weyl Curvature 

The gauge freedom can makes the congruence of n a -integral curves shrink into a point or take away to 
infinity in the direction of n a . Indeed, when we take an affine gauge the null generators are complete 
in affine parameterization by definition and go away to infinity. Nevertheless, since we want to study 
the geometrical structure around + in n°-limit, it is essential to choose a compact gauge il' = u>Q, in 
which the n°-integral curves are incomplete in affine parameterization (possible to compactify) though is 
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complete in a original parameterization n a — (d/du') a . Under this gauge condition, we discuss the limit 
structure of the congruence of the rt a -integral curves in n a -limit. 

When we choose the compact gauge we should be careful about the angular dependence of u, since 
we include the case of vanishing lo. If we allow lo to have angular-depending irregularity in n a -limit, 
the irregular angular dependence may change the topological structure of the congruence of n a -curve. 
Because a conformal transformation should be diffeomorphism, we must choose such lo that is smooth 
with regard to the angular and positive. 

For example, let us consider a sphere. The metric on a sphere with unit radius is ds 2 = d9 2 + sin 2 9d<p 2 . 
When we perform the conformal transformation such that a conformal factor vanishes at 9 = n/2 (for 
example, lo = cos 9, its metric is ds 2 — cos 2 9(d0 2 + sin 2 9d<p 2 ) and is not diffeomorphic to the original 
sphere. We must choose a regular gauge to avoid such a case. Then we never allow that dependence for 
the gauge transformation f2' = loQ from an original affine gauge. 

Newman-Penrose equation^HI relates the shear of the congruence around to Weyl curvature as 
the following; 

-n a V a A + m a V Q ^ = + (37-7)A 

+ (3a + /3 + 7r-f)j/ + ^ (4) . (28) 

We expand each term of the above equation in powers of and leave the only leading term. From 
(fTHfl the order of A is 

A = 0(0). (29) 

The spin connection v can be written as 

v= -im a V a fc (1) + O(0 2 ) (30) 

by the equations (0) and J7J). We can choose satisfying 

m a V Q fc (1) = (31) 

from the relation 1)25(1. For we can write m a as m a = m e (d/d8) a + im^(d/d(f)) a using two spacelike 
vectors, (d/d9) a , (d/d(j>) a orthogonal to n a ,l a . Then if we transform to from affine gauge into another 
gauge, from the equation l|24(l . the equation l(31|l becomes 

d f n a V a u \ _ d f n a f a oj \_ n 

We may choose ui so that n°V a 0/0 2 does not depend on 9, (f> and is positive everywhere. Therefore, we 
see that for satisfying 1(31(1 . the order of v is 

v = O(0 2 ). (33) 

Finally, let us compute the order of the spin connection 7 defined by 7 = l/2(n a n b V a h — n a fh h V a l{,). 
The real part of 7 is 

7 + 7 = n a n b VJ b 



m + 0(0) (34) 
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from the equations and Q14JI. Since the imaginary part of 7 becomes 

7 -7 = -n a m b V a m b , (35) 

it depends on the direction of m a which we have not yet determined. Now let us determine the direction 
of m a as follows. We determine the direction of m a so that on it will be parallelly transported along 
the null geodesic generators on J? + . That is, 

n a V a m b = on J?+. (36) 

Moreover in the direction away off J^ + , it is parallelly transported along —l a , that is, 

rV a m b = everywhere. (37) 

Thus, 7 — 7 becomes 

7-7 = 0(0). (38) 
So we see that from the equations l|34|) and (|38|l . the order of 7 is 

7= M +O (0). (39) 
From (|29[l . (|33|) and l|39|l . we can write the equation (|28|) as 

n a V a X {1) = 3/i (0 )A ( i) - * 4 (i) (40) 

in the leading order. 

Integrating this, we obtain the relation 



(41) 



A (1) = - exp ^3 J ^( )(u")du"j • J du" * 4(1) (u") exp ^-3 J fj.^(u'")du 
where u' is the parameter of n' a and is defined by n' a = (d/du') a . 

4 Topology of Event Horizon 

Now we discuss the relation between the topology of event horizon and the n a -limit structure of n a - 
integral curve congruence. Pseudo-Newman-Penrose basis provides the following timelike submanifold 
near ^ + but not on In this region, since u> is monotonically decreasing function in the direction 

of n a supposing sufficiently large parameters, n a is timelike from ij9")l l|24|) (|25|) and orthogonal to a spatial 
hypersurface C e which is defined by small constant e as f! = e < 1. From the smoothness of the 
pseudo-Newman- Penrose basis, a sufficiently large sphere S e is defined on C e and spanned by m° and rh a . 

n' a = (d/du') a generates a one-parameter family of diffeomorphisms ip u i on M, and we have a family 
of two dimensional surface S e (u') = ^v(«S e ) and a submanifold 

We = {p£tp u '{S e )\u' £ (-00,00)}, (42) 

where "00 identically maps S e onto itself. 

Since M is not a manifold at i + , we want to discuss the structure of S e (u') without dealing with z + , 
that is within openset Nevertheless, the inextendible timelike curves on W e does not have its 
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Figure 1: C e is the spacelike surface where fi is constant. On that a sufficiently large sphere S e is given. 
n a generates a timelike submanifold W e from it. It determine a two boundaries of its past set I~{W e )- 
Inner one is Hi and outer one is H . 

endpoint on the future null infinity or rather on the outside of For sufficiently small 

e and large S e , there is W e outside of black hole, since the black hole region M\I~ (J? + ) is compact and 
n a n a = n a V a f2 = fcmfi + ... < means decreasing of in n a direction. 

Since <S e (0) = S e is topologically a sphere, I~(J? + ) is decomposed by W e into an inside subset A and 
an outside subset B (W e C A). If (B n M) n n A/) 7^ 0, W e can cross the event horizon H+ . 

Then 

(BnM)n(r(/ + )nM) = (43) 
r(w e ) = r(/+)nM (44) 

The boundary of the causal past of W e is composed of ingoing part Hi — A\ W e generated by ingoing 
past null geodesies and outgoing part H = B \ W e generated by outgoing past null geodesies. From the 
continuity of M, the ingoing part approaches to the event horizon. Considering liniu'—Kx, S e (u'), we may 
observe the event horizon in n a -limit H + n Ue is possibly not homeomorphic to the cross section of J^ + , 
where Ue is neighborhood of E in M. 

On the submanifold W e we will observe the deformation of the section S e {v!) as a deviation of the 
geodesic congruence of n a . Especially the index theorem of the Morse theory [Tj^ is powerful to see 
topological aspects. Next, we mention the method to examine the limit of S e near J? + . The following 
corollary |lfi| of the Poincare-Hopf 's theorem is useful. 

Theorem . Let M be a compact n- dimensional (n > 2 is an odd number) C r (r > 1) manifold with 
Si U £2 = M and Si n £2 = 0- X is any C r ~ 1 vector field with at most a finite number of zeros, 
satisfying the following two conditions: (a)The zeros of X are contained in Int M. (b)X has inward 
directions at Si and outward directions at S2. Then the sum of the indices of X at all its zeros is related 
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to the Euler numbers o/Si and £2' 



X (S 2 ) - = 2 indexpO, (45) 

where index(JT) is given by the alternating sum of the Morse number fi^ as index = X)k( — ^) k / i k- The 
Morse number fj,k is the number of a critical point (zero of vector field X) whose index is k. The index of 
a critical point is given by the number of negative eigenvalue of Hesse matrix H a b — V a X{>. This theorem 
means that when the tnnnlocv of slices nf a. manifold chances, there must he a. zero of the vector field on 
it and eqi ^] 

Here \ ic. 
As discus Q : is 

an embed 

16) 

in which . 




Figure 2: Timelike submanifold W e is diffeomorphically embedded into a trouser-like manifold Af. On 
the boundary of W e there are zeros of the tangent n la . 

Because W e and Af are a three-dimensional manifold, we apply the theorem to Af. Let us consider a 
continuous tangent vector field X of Af such that X is the natural extension of <fi*(n a ) 

X% {We) =<p*(n a ), (47) 
is outward directed on upper boundary S 2 and do not have a zero on Af \ (j)(W e ), since Af\ 0(W e ) = 

S 2 X (-C3O,0]. 

Here, we should pay attention to the fact that, in this case, zeros can exist only on the boundary 
of W e , since W e is complete about its own tangent n a by definition. As easily understood, the upper 
boundary of </>(W e ), 3(0(W e ))\ Si is the set of endpoints of X a since n a is complete on W e . Nevertheless, 
the endpoints (zeros) with index= do not affect the theorem and are extendible endpoints since Morse's 
lemma [19] states that the Morse critical point is generic. By small reparametrization, we have X' a 



11 



which is almost same as X a except for small neighborhood of d(4>(W e )) \ Si, and have only one zero on 
d(<j){We)) \ Si with index= 1. 

Therefore we have to be careful to apply the theorem. Nevertheless, it would be possible to relate the 
index of the theorem and the congruence of null generators around the zero. From Morse's lemma, we 
suppose that the zero of the null generators is isolated. Since M is three dimensions, the behavior of the 
congruence around the zero which is on W e is illustrated in FIG El 

The left figure in FIG|31is a three dimensional conformal embedding with one dimension suppressed. 
The right figure in FIG 01 represents the close-up figure of encircled points of the left figure and makes 
clear the behavior in the vector field X a and n a in the neighborhood of their zero. We should note that 
though in the left figure in FIGOHtV is drawn as two dimensional surface rather than three dimensional 
surface, in the right figure the neighborhood of the zero is drawn as the neighborhood of the zero on 
three dimensional surface. In this figure, there is a zero with index = +1. In the neighborhood of the 
zero, we can introduce a coordinate system (xi, X2, X3) on J\f. A heavy arrow and a dashed heavy arrow 
are on xi axes, dashed arrows are on a^a^-planes, dotted arrows lie on xia^-planes and solid arrows lie 
on iriiz^-planes. In the left figure of FIGEl the heavy arrow is on the front side while the dashed heavy 
arrow is on its reverse side. 

Seeing the region contained in 0(yV e ), X2, x^)\ — x\ +a;|+cc§ > 0}, we can extend its neighborhood 
into the other region in M \ d{<j)(W e )) 1 {(x%, X2, %3)\ — x\ + x\ — x\ < 0}. For it contains the region 
{(x%, X2, £3)1 — x\ + x\ + x\ > 0} but does not contain the other region {(x\, X2, x^) \ — x\ + x\ + x\ < 0} 
in a certain coordinate system. The limit behavior of the shear and the expansion of the vector field 
n a = g ab V b Vl tells us information about its zero since V b n a = Vf,V a il gives Hesse matrix. Approaching 
to that zero along the heavy arrow, the dotted arrows on the 2;ia;2-plane in FIGEHare going away from the 
zero, and dashed arrows on the xia^-plane in FIG|21are approaching to the zero. Therefore, approaching 
to the zero along the heavy arrow in FIG|21 the shear of the vector field become larger than the expansion. 

Note that because at the zero with index= 1 tangent vectors X a vanishes, it is Lifshitz continuous 
(that is, at least C 1_ ) at that point. Then, the embedding of W e into Af is C . We should note also that 
when we consider the embedding of a neighborhood of that zero, the into-mapping is not self-intersecting 
because of the fact we introduce a trivial patch with R 3 in the neighborhood of that zero. Furthermore 
we comment that because in Schwarzshild spacetime, there are nothing but zeros with index= of vector 
field X a at the boundary of W e . W e have a trivial completion into JV ~ S 2 x [0, 1]. 

To summarize, for the Hesse matrix H a b = V a nb = V a Vf,r2, we solve the eigenvalue equation 

H ab x a = Ax b = Ah ab x a (48) 



where h ab is the metric of W e . Considering three basis vectors of W e , n a = n a / y/\n c n c \, (m a ± m a )/\/2, 
the eigenvalues of H ab are given by 

n a n b V a n b 

/ri — ] , M±|A| ■ (49) 

In compact gauge, n' c n' c is negative not on ^ + and converges to zero at v! — » 00. Therefore the first 
eigenvalue is negative and the index is determined by inequality < |A|. In this gauge, SI' = uiQ, is 
decreasing positive function about v! and converges to zero in u 1 — > 00. Then we are able to choose two 
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gauge conditions, they are d u 'Uj/uj — ► and d u iui/u) — > /3 < and corresponds to w oc (a > 0) and 

w oc e"' 3 " . lu oc e~^ u , however, means n a n a is not zero in v! — > oo and this is a bad aspent. When we 
take this gauge, the n a -integral curves go into a timelikc infinity. Nevertheless the timelike infinitiy is 
generally pathological and M is not manifold there. It will mislead topological investigations. 
Therefore we must choose a gauge 



a 



(50) 



where higher order of // about f2 is eliminated by higher order gauge condition of w. 

So if the behavior of Weyl curvature on J? + is ^4(1) ~ u'~ e as u' — > 00, A' behaves as A' oc u' l ~ e Q! 
by eq. (|41|) . The ratio of the shear to the expansion is 

A' 



From 



e 0' + O(0' 2 ), 



(51) 



> n' c n' = 



1 v 1 



oc u 



n ' c w c n' = fc (1) 0' + 

/-2a 



(52) 
(53) 
(54) 



where it shoul 




Figure 3: The behavior of the vector field n a = g ab Wb^l around the zero on a trouser-like manifold. The 
left figure is Af with one dimension suppressed and the zero is encircled. The right figure is close-up 
of the encircled point of the left figure. Each heavy arrow is in the direction of X2-axis. Each dashed 
arrow lies on the a^a^-plane and each dotted arrow lies on the sia^-plane. We should note that in 
the left figure Af is drawn as two dimensional surface (actually it is a three dimensional surface) but 
that each right figure is as drawn the neighborhood of a zero on three dimensional surface. The region 
{ (xi , X2 , £3 1 — x\ + x\ + x\ < 0)} is not contained in W e and is the small extension in the direction of 
X a . When we consider the expansion and the shear of the vector field n a along heavy arrows (the heavy 
arrow and the dotted heavy arrow), then we observe that approaching to the zero along heavy arrows, 
the shear becomes larger than the expansion. 
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that the 



Therefore, if there is a satisfying 2 — 2a — e > 0, this ratio can diverge, that is, the shear can be much 
larger than the expansion in this gauge. As shown above, though even the topological feature changed 
by the gauge condition it is guaranteed by the existence of a gauge providing not decaying 
event horizon is topologically non-trivial like FIG|5] 

The point is that if Weyl curvature falls off later than some power of parameter of tangent n a of J? + , 
two black holes will not coalesce eternally. Now that we obtain the sufficient preparation, we can restate 
our sufficient condition for not coalescing black holes as follows; 

Condition 1 (Sufficient condition for separated two black holes). When we transform the con- 
formal factor from affine gauge into another gauge — > 0' = 00, the number of black holes is not only 
one, if there exists an unphysical spacetime (M, 2 g a b) with the conformal factor satisfying the below 
conditions. 

(1) There is a real number a such that = 0(9, (j))u'~ a and liniu^oo = 0, where u and u' are parameters 
of the null geodesies on J? + in affine gauge and in another gauge, respectively and 0(6.<f>) is the smooth 
function on a sphere which is positive and not singular everywhere. 

(2) There is a real number a > satisfying the following condition; For any number L > 0, there are 
positive numbers K such that if u' > K, then \u' 2 ~ 2a ~ t \ > L, where e is an exponent appearing in ^4(1) 
extended in the power of u' , that is, ^4(1) ~ u ' 



l—e 



Finally also we note that Einstein equation is only used through the definition of asymptotic flatness 
in the section 2 and Newman-Penrose equations. The condition 7 in Definition ^ m the section 2 is equal 
to the condition that energy momentum tensor goes to zero at future null infinity rapidly if we impose 
the condition that Einstein equation holds. We used this condition in the ea. l|12(l and so on. (Though 
in the book of Wald (reference |17p vacuum Einstein equation is assumed in the neighborhood of future 
null infinity, we weaken this requirement.) We never suppose something like energy conditions. 

Since this condition is a sufficient condition, it is not clear whether the spacetime satisfying the 
condition is realistic. In the next section, by proving in C-mctric the <|54ll diverges (e < 2) in appropriate 
gauge, we show that the event horizon is not that of a single black hole but that of two black hole 
connected at the virtual boundary of conformally compactified manifold. This is what we already know 
about the causal structure of C-metric spacetime. 

5 C-metric as an concrete example 

To verify the relevance of the condition, we take an example which satisfies the above condition. Now 
we treat the vacuum C-metric [21], [22 as an example of the spacetime in which two black holes never 
coalesce from the viewpoint of an observer outside black holes |20| . The line element of the vacuum 
C-metric is written in the form of 

ds 2 = -Hdu 2 + H~ 1 dr 2 + 2Ar 2 H~ 1 drdx 

+ (G~ 1 + A 2 r 2 H- 1 )r 2 dx 2 + r 2 Gdip 2 , (55) 
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where 



H = 1 - 2mr~ 1 + 6Amx + ArG, x - A 2 r 2 G, (56) 
G = 1 - x 2 - 2Amx 3 . (57) 

In BMS coordinate[2B], it is expanded at large distance. Since in our discussion it is sufficient to consider 
the only neighborhood of J^ + , we may transform the metric <|55|) to a BMS coordinate written by 



where four functions are 



ds 2 = -fdu 2 - 2e 2f3 dudr - 2Udud9 

+r 2 (e 2 "<d9 2 + e~ 2 ~< sin 2 9dcj) 2 ), (58) 



f(u,r, 9) = l-« + ofi), (59) 

(3(u,r,6) = ~%^ + 0(l), (60) 

7 (u,r,9) = ±j3. + (l), (61) 

U(u,r,0) = -(c e + 2ccot9) + 0[ - ), (62) 



and M(u, 9) is given by 



M u = -c 2 u + -(c ee + 3c g cot 9- 2c),„. (63) 

We should note that these four functions are expressed by one function c(u, 9) whose differentiation by it 
is called the news function, 

— 3a 3 m sin 2 9 

c - = IZr^r- (64) 

We introduce a new coordinate by 

v = u + 2r. (65) 

In this coordinate, the metric Ij58(l becomes the following form, 

ds 2 = -(/ - e 2f3 )du 2 - e 2f) dudv -2Udud9 

+-(v - u) 2 {e 2 " / d9 2 + e~ 27 sin 2 9d<p 2 ). (66) 

Furthermore, we introduce a new coordinate V — l/v so that infinity along outgoing null geodesies will 
correspond to V — 0. The metric components in the new coordinate (u, V, 9, <p) are given by 

,2/3 



ds 2 = -if - e 2f} )du 2 + %-dudV -2Udud9 

1 



+ -(v - u) 2 {e 2 ^d9 2 + e- 2 ~< sin 2 d(f> 2 ). (67) 



In turn, we conformally transform the above physical metric into 

ds 2 = -(/ -e 2p )0 2 du 2 + e 20 dud0-2U0 2 dud9 



+-(1 - 0u) 2 {e 2 ~<d9 2 + e- 2 ~< sin 2 9dcf> 2 ), (68) 
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where this unphysical metric is related to the physical metric by g a b — 2 g a b and = V. The four 
functions (|5*9")) . (|5T7|) . and are expanded as 

/ = 1 -4M(u,6>)0 + O(0 2 ), (69) 
13 = -c 2 (u,9)0 2 + O(0 3 ), (70) 
7 = 2c(u,6)0 + O(0 2 ), (71) 
J7 = -(c g + 2c cot 6<) + [-2c(c<? + 2c cot 9) 

+2{2c u + 3cc e + Ac 4 cot 6)}0 + O{0 2 ) (72) 

in powers of conformal factor 0. 

Now let us determine the basis vectors defined in section [3] from the metric (|68f) . It is easy to 
decompose the metric <|68[1 into the form of 

9ab = -(fyi)0 H )Wi> - 2/ (a ?i b) + 2m (a m &) . (73) 

We see that the four dual vectors, 

n a = (d0) O) (74) 
/a - -ie 2/3 (^)a, (75) 



1 / 4C/0 2 



2V2 y V (i-w0) 2 

+e 7 (d6>) a -ze _7 sin6>(d(/>) a ), (76) 



1 /-, ^/ --v 4C/ 2 ,J X 
+e 7 (d6») a + ie ~ 7 sin ^c^a), (77) 



satisfies the eauation (|3T)l . 

We raise the indices of the dual vectors to obtain four basis vectors, 

n a = 2e-w(^y 
\ ou . 



+ [ 4(/-e^)e-^ + 16e-^ IT ^ F ]( <V V 



(1 -u0) 2J V30 

(81) 



l-u0\dO) sin 6»(1 - U 0) V d<j> ) ' 
V2e^ f d\ a . V2e 7 fd\ a 



l-u0\de) sin0(l-u0) 
From (|74|l and l|78(l . we compute the norm of n a as follows. 



n a n a = 4(/ - e^e^ 2 + lSe- 4 *- 2 ^^ 20 * . (82) 

(1 — U0) Z 
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Since the first order of f2 in the equation (|82[1 vanishes, we see that the choice of conformal factor = V 
is afhne gauge at J? + . That is, in the choice of this gauge 

M(o) = ~P(o) = _ 2 fc W = ^ 83 ^ 

holds. The shear A in this gauge is 

A = —m a rh b V a nb 

= -4c»0 + O(0 2 ) (84) 
3a 3 TOsin 2 ^^. ,„_. 

where we used the fact the news function of the vacuum C-metric is given by ljfi4|l and u is the affinc 
parameter of the null generator of J^ + . 

For the investigation in previous sections we require a gauge transformation from the afhne gauge into 
the gauge such that ,J? + can be compactified in future direction. Then we examine the ratio of the shear 
to the expansion of the congruence of n a -integral curves which plunge into the zero at n°-limit. In affine 
gauge, the shear in this direction behaves 

A~ A (1) r}~ J_ ' ( 86 ) 
u 

as u — > oo and £1 — * 0. Under the transformation 0^0' = 00, there is the relation between affinc 
parameter and the parameter of the geodesic in new gauge as follows 

u 'a+l 

u = —— (a^-1) (87) 
a + 1 

where we set = u'~ a (a > 0). 

From J23 {SU (JHnj and 1 t ne ra ti° 01 tne shear to the expansion is given by 

— - u' a+1 & - u n - a . (88) 

Then, we see that when u — > oo{v! — > oo), this ratio diverges to infinity for < a < 1. We ensured 
the existence of the conformal factor 0' = 00 such that the n a -limit ratio of the shear to the expansion 
of the n a -congruence around diverges. We also find easily this gauge is compact gauge, 

lim = 0. (89) 

u — >oo 

On the other hand, though Newman-Penrose formalism is ,of course, true, we confirm that the behavior 
of Weyl curvature is obtained directly from the metric 

1 



*4(i) - (90) 

If we choose lj so that < a < 1, 0<e<2 holds. Therefore, two conditions (1) and (2) of the sufficient 
condition stated in the end of section 4 are satisfied. This implies that the two black holes of the C-metric 
are eternally separated. 
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6 Summary and Discussions 



We have investigated the limit structure of the congruence of the n a -integral curves and its orthogonal 
spatial section introducing a pseudo-Newman Penrose basis vectors. It was shown the sufficient condition 
in order that two black holes never coalesce to an asymptotic observer outside black holes in such a 
spacetime with two black holes. This is that if the tail of outgoing gravitational wave ^4 remains at 
sufficiently late time for an asymptotic observer, two black holes remain separated forever. 

In other words, we clarify that the behavior of Weyl curvature ^4 in the neighborhood of future null 
infinity at late time rules the geometrical structure of the congruence of tangent of J^ + . Especially it is 
considered in spacetimes with two black holes, the two black holes separated forever is caused by the shear 
of the congruence through the Poincare-Hopf theorem. Moreover, we made sure that Weyl curvature of 
C-metric remained so late that it satisfy the sufficient condition. 

Here we comment that to discuss the final stage of the gravitational collapse we have only consider 
causal structures in the limit of the future direction instead of supposing stationary spacetime. Though 
at first sight they might correspond to each other, the authors know nothing about that. More advanced 
studies in the course of the present investigation will tell something about that. 

Studying black hole perturbation in some exact solution background, it have been known that there is 
gravitational wave tail in late time limit [?]. For example, it is easily confirmed that outgoing gravitational 
wave is far smaller than the late discussed in this article for Schwarzschild spacetime. This is valid since 
Schwarzschild spacetime is a single black hole spacetime. We guess that also in other single black hole 
spacetimes the situation is same. Otherwise, the topology of black hole will be changed by the effect of 
the gravitational wave perturbation. 

As a future work we need investigate the physical and geometrical meaning of the parameter vl ', 
althogh in this article, we treat it as the parameter such as future null infinity is compactified. It is 
interesting to investigate the relation between the compact gauge and the BMS group, however, at this 
stage, we do not know the answer. 

While we concentrate on the gravitational radiation "J 4 of Weyl curvature as leading order contribu- 
tion, there are possibilities that other component of curvature may contribute to the strucutre of tangent 
n a of .y + in higher order. For example, in Majumdar-Papapetrou spacetime, while gravitational field 
radiation $4 vanishes, we observe that ^2 contributes to the congruences of n a in a higher order. 

Furthermore, when we consider a higher dimensional black hole spacetime, a investigation like in the 
present article will reveal the mechanism of topological black holes, since a higher dimensional counterpart 
of the C-metric spacetime is black ring solution with vanishing angular momentum|ll| . 

Using the result of this article, what can we know about the collapsing stars from the observation 
of the gravitational radiation ? The sufficient condition for separated black holes suggests the power of 
late time tail radiation will let us know about the final number of black holes. Of course, to get powerful 
statement about that, we should make clear the astrophysical meaning of the gauge choice used in this 
article. 



18 



Acknowledgments 



We thank Dr. Nakamura for useful comments. We would like to thank Professor A. Hosoya and T. 
Shiromizu for continuous encouragement. This work was supported by a 21st Century COE Program at 
Tokyo Tech "Nanometer-Scale Quantum Physics" by the Ministry of Education, Culture, Sports, Science 
and Technology. 



References 

[1] W. Israel, Phys. Rev. D164, 1776 (1967) 

[2] W. Israel, Commun. Math. Phys. 8, 245 (1968) 

[3] B. Carter, Phys. Rev. Lett. 26, 331 (1971) 

[4] D. C. Robinson, Phys. Rev. D10, 458 (1974) 

[5] D. C. Robinson, Phys. Rev. Lett. 34, 905 (1975) 

[6] S. W. Hawking, Commun. Math.Phys. 25, 152 (1972) 

[7] P. T. Chrusciel and R. M. Wald, Class. Quant. Grav. 11(1994) L147 

[8] D. Gannon, Gen. Relativ. Gravit. 7 (1976) 219, J. L. Fricdmann, K. Schleich and D. M. Witt, Phys. 
Rev. Lett. 71 (1993) 1486, T. Jacobson and S. Venkataramani, Class. Quantum Grav. 12 (1995) 
1055, S. Browdy and G. J. Galloway, J. Math. Phys. 36 (1995) 4952, G.J. Galloway, K. Schleich, 
D.M. Witt, E. Woolgar, Phys.Rev. D60 104039, (1999) 

[9] S. A. Hughes, C. R. Keeton, P. Walker, K. Walsh, S. L. Shapiro and S. A. Teukolsky, Phys. Rev. 
D49 (1994) 4004, A. M. Abrahams, G. B. Cook, S. L. Shapiro and S. A. Teukolsky Phys. Rev. 
D49 (1994) 5153, S. L. Shapiro, S. A. Teukolsky and J. Winicour Phys. Rev. D52 (1995) 6982, P. 
Anninos, D. Bernstein, S, Brandt, J. Libson, J. Masso, E. Seidel, L. Smarr, W. Suen, and P. Walker, 
Phys. Rev. Lett. 74 (1995) 630, S. Husa and J. Wmicour,Ptys.Rev.D60 084019 (1999) 

[10] M. Siino, Phys.Rev. D58, 104016 (1998),M. Siino, Phys.Rev. D59, 064006 (1999) 

[11] R. Emparan, H. S. Reall, Phys. Rev. Lett. 88, 101101 (2002) 

[12] R. Penrose, Proc. Roy. Soc. Lond.,A284, 159 (1965), R. P. Geroch, "Asymptotic structure of space- 
time" in Asymptotic structure of space-time,ed. F.P.Esposito and LWitten (New York:Plemium). 

[13] A. Tomimatsu, Phys. Rev. D57, 2613(1998) 

[14] J. B. Hartlc and S. W. Hawking, Commun. Math. Phys. 26, 87 (1972) 



[15] V. Pravda and A. Pravdoda, gr-qc/0003067 



[16] R. D. Sorkin, Phys. Rev. D33.978(1985) 

[17] R. M. Wald, General Relativity University of Chicago Press, Chicago, 1984. 



19 



[18] E.Newman and R.Penrose J.Math.Phys. 3, 566 (1962) 

[19] J. Milnor Morse Tiieory,Princeton University Press, 1963 

[20] F. H. J. Cornish and W. J. Uttley, Gen. Rel. Grav. 27,439(1995) 

[21] W. B. Bonnor, Gen. Rel. Grav. 15,535(1981) 

[22] W. B. Bonnor, Gen. Rel. Grav. 16,269(1983) 

[23] H. Bondi, M. G. J. van der Burg and A. W. K. Metzner, Proc. R. Soc. London A269,21 (1962) 

[24] P. O. Mazur, J. Phys. A15, 3173, (1982), Acta Phys. Polon. B14, 219 (1983) 



20 



